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ATMOSPHERIC DENSITY MODELS 
by 

Alan C. Mueller 


1 . 0 Inc roduct ion 

To compute the density at any point in the atmosphere, 
the Mission Control Center orbit prediction proRram will 
employ a model developed by Jacchia (Reference 1). Although 
the model is quite accurate, it requires a large amount of 
computer storage and execution time. Because of software 
constraints, these undesirable features make the Jacchia 
model ill-suited for the Shuttle onboard program. However 
if a different model is to be used on board, mission con- 
straints require a c>irtain degree of compatibility with the 
Jacchia model used on the ground. Thus the intent of this 
study is to develop a simple atmo-ipheric density model that 
simulates the Jacchia model . 


2.0 Atmospheric Density Models 

To date, there exist only two analytical, dyn.'imic models 
which include most or all of the important characteristics of 
the earth's atmospheric density. One is the Jacchia model 
mentioned in the introduction and the other is a model devel- 
oped in the U.S.S.R. (Reference 2), The USSR model was used 
during ASTP , the joint Russian-Atnerican mission. The follow- 
ing characteristics of the density are included in both models 

1. Variation with solar activity (dally and 
3 month average) 

2. Diurnal variation 

3. Variation with geomagnetic activity (yearly 
mean and 3 hourly index) 
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4. Semi-annual variation 

5. Variation with height 

In the Jacchia model, the first four variations are com- 
puted as char.^-^s in the e?:ospheric temperature Too • One may 
then obtain a rather complicated solution to the diffuse equi- 
librium equation, knowing Too and the height h . The density 
P may be determined from the solution of this equation for 
the different gas constituents. In addition the Jat.chia model 
also computes the variation in density due to a seasonal - 
latitudinal effect which becomes important for altitudes below 
300 Km. 

The USSR model, on the other hand, computes all the above 
characteristics directly as variations in the density and is 
computationally much simpler that the Jacchia model. The USSR 
model does not account for the seasonal - latitudinal vari- 
ations. Unlike the Jacchia model, the coefficients used in the 
USSR model are not constants but are a function of the chang- 
ing solar activity level. The coefficients are given In a tabu- 
lar form for different values of the 10.7 cm solar flux 
(Fio j) .a measure of the solar activity. Unfortunately, the 
coefficients are not known for the expected values of ^ 

in the early 1980' s which renders the model useless for many 
shuttle flights. Yet tests have shown that the USSR model 
agrees with the Jacchia model for solar activity levels V 
which the coefficients are known (Section 2.2.2, Reference 3). 
Therefore the USSR model is certainly a candidate for the on- 
board program if the coefficients can be determined for the 
expected values of F^^ ^ during the Shuttle flights. 

One approach to determine the coefficients would be by a 
least square data fit using some iceration scheme such as the 
Davidon technique to adjust the coefficients. The data to be 
fit could be obtained from the Jacchia model. Although this 
type of approach could be extremely accurate it also would 
require a considerable theoretical and computational effort. 
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AMDB, u density model developed by Gus Babb (Section 
2.2.4, Reference 3) may offer some insight In the problem of 
determininR coefficients for the USSR model. The AMDB model 
does account for the height and diurnal variations of the 
density model. The other variations are included implicitly 
by calibration to the Jacchia model. As in the USSR model, 
the coefficients in AMDB, are implicit functions of the solar 
activity level, but unlike the USSR model, AMDB coefficients 
may be determined easily by a procedure known as calibration. 
The A.MDB model's simple formulation allows one to Invert the 
the problem; given the densities at different points in the 
atmosphere one can explicitly solve for the coefficients. 
However, the AMDB model does net explicitly account for all 
the variations and has been found to be not as accurate as 
the USSR model (Section 2.2.5, Reference 3) 

Thus a new approach is proposed which is to develop a new 
density model using the accuracy advantages of the USSR model 
and vhe calibration advantages of AMDB. With this new di- 
rf^ction, it may be possible to reach the goals established in 
the introduction. 

3 .0 Comparison of A.MDB and USSR models 

Before attempting to develop a hybrid of the USSR and 
AMDB mf)dels it is best to understand the similarities and 
differences between the two models. Since the ground Jacchia 
model will neglect the short period variations in the solar 
activity level and geomagnetic index (because of their small 
effect ar.d for operational reasons), these variations will 
not be included in the comparison. With these effects ne- 
glected, the USSR model reduces to the following form for the 
density p 
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where p is the ni^ht time vertical profile given by 
n 

- exp (a| - a^ - a^' ) , (21 


Kj is the variation duo to the diurnal effect given by 


K, » 1 + B(h)-b (a ,6 ,a ,6 ) 
1 V V 8 8 


(3) 


and 

by 


is the variation due to the semi-annual effect given 


K 2 - 1 + D(h).E(d) 


(4) 


The function B is a rather complicated function of the alti- 
tude above the oblate earth h . D is a linear function of 
the height and b is a trignometric function of the right 
ascension and declination of the vehicle (a ,6 ) and that 

V V 

of the sun (a ,6 ) . And finally E(d) is a tabular 

8 S 

function of the day of the year. 

It Is the diurnal term (equation (3)) which prevents the 
USSR model from being analytically inverted (given the density 
solve for the coefficients in the model). Of the 15 coeffi- 
cients used in this model, 10 are contained in the complicated 
expiession for tlie diurnal term. 

In the AMDB model the density is given by 

p = p^ exp (F + G) (5) 


where F is the dusk density vertical profile given by 


F -■ Uj + a^h + a^/h 


( 6 ) 


and G is the diurnal term 

G = (bi + b2h + b^/h)-g ' “s ’ ^ 


(7) 
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whoro Uj, , a^, b^, , b^ are model parameters. 

The major dil'lerence between the nxidels is the fact that 
the diurnal term it» AMDB is an additive term within the ex- 
ponent. This results in an adequate but simple form for the 
maRnitude of the diurnal term as a function of the heiRht, 

The trignometric function c , however, does not seem to accu- 
rately model the angular dependence of the diurnal term. 


4.0 The Diurnal Term 

From the comparison of the two models it is apparent that 
the diurnal term results in complications in both models. In 
the AMDB model it is inaccurate and USSR model it is too com- 
pl j cated . 

Assume that the density at a specific height but at 
arbitra-y angular coordinates is given by 


P = (1+Bg*) 


0 < g* < 1 

B > 0 


( 8 ) 


such that when g* = O , the density is at a minimum and 
R* = 1 the density is at a maximum. B then gives the dif- 
ference of the densities at two points at the same height: 
one directly in the atmospheric; bulge; where the density is 
greatest and the other point on the opposite end of the line 
which passes through the bulge and the center of the earth 
where the density is lc>ast. Assume also that the function 
g* varies by the following 




L m 
' D ^ 


(9) 


where L is the distance between the vehicle and the least 
dense point, D is the distance between the least and most 
dense points and m will be defined later. Sinct^ these points 
are restricted to a fixed altitude. D is equal to two times 
the satellite's distance from the center of the earth. I) = 2R . 
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The distance L is given by 

L - /(x-x + (y-v + (z-z )^' 

in m ID 


where 


x,y,z define the vehicle position 

X , y , z define the least dense point 
m m m 


X 

- R 

aoB 

( Y + ‘") 

008 ( -6 ) 

s 

- R 

aoB 

y ooB t 

m 





8 




B 


- R 

sin 

( Y + ti ) 

aoa { 


s 

- R 

sin 

y OOB 6 

m 





S 




s 

z 

= R 

sin 

(-6 ) 

* - R sin 

6 





ms 3 


( 10 ) 


( 11 ) 


Y * a + t 
s 


( 12 ) 


: is the lag of the bulge behind the sun because of 

rotation of the eartn. 

Inserting all this into equation (10) and simplifying, one 
finds 


L = 2R ( 


1 -*• coa ij; .4" 


where 


aoe li' = p 


z 3 in 6 + cos 6 (x • ’c y + y sin y) 

s s 


] 


(13) 


and thus 


m 


m 


«* = ( g ) = ( 


l+OOS'i) ,2 
o ) 


(14) 


Finally by trignomotric identities onf' obtains 


g* = IJ 


(15) 
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By careful obKorvation of the diurnal function In the USSR 
one finds it has the form 


b 


a03 + c 708 — 


2 


2 


( 16 ) 


where is the angle formed by the least dense point and 

the vehicle point and the center of the earth. 
angle for the most dense point. But <• is always small 
(|c| < .05) and if this term is neglected the angular func- 
tion in the diurnal term is exactly that obtained in equation 
(15). With all this in mind we now attempt to define a hybrid 
AMDB/USSR model. 


5.0 Development of Hybrid Model 

From considerations ot sections 3 and 4 the following 
model is proposed 


p* = exp ( A + B ) 


where A is the night time vertical profile 


A 


+ 


a^h + 



and B is the diurnal effect 

„ . , b, * b^h * ^ 


m 

l+apaj) V n 
2 ^ * 


(17) 


( 18) 


(19) 


where aoe 4 < is given by equation (13). 

The seasonal-latitudinal variation may be included in 
this new model. Since the Jacchia model computes this term 
as an explicit variation in the density one may incorporate 
this in the new model easily. The seasonal -latitudi nal vari- 
ation that is given by .Jacchia is 


= n ♦ 


10 ' 
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who ro 


c 


6 

V 

(0.02) (h-90) — exp 


• oin 


360 

(d+100) 

y 


-0.045 (h-90) j 6^ 

( 21 ) 


/ 


h is (riven in km, d Is the number of days into the year 
and Y is the number of days in a year. This term may bo 
rewritten so that the final form ol the model Is 


p ■ p exp (A+B+C) 
o 


( 22 ) 


where 


C = .01605 ♦ I z j • — (h-90) (?xp 


-0.045 (h-90) 




• oin 


360 

(d+100) 

Y 


(23) 


The term underlined may be computed once with the initial day 
d and days in the year Y and then assumed constant . 

6 .0 Calibration Procedure 

If one assumes that the coefficient m in equation (14), 

p and the bulge angle 4> are known constants and unaffectea 
o 

by the solar activity or geomagnetic changes then the hybrid 
model has 6 coefficients which are determined through calibra- 
tion to the Jacchia model. 


Let us first Rive the K^fneral solution to 3 etjuations 
in 3 unknowns with the following form 


Cj * Cj * C3/W, . f| 


i » 1 . 2 , 3 , ( 24 ) 


The solution of these equations is 


W 1 .W 3 




(Wj-W^) 


Wi W3 


w,W2 


^2*^1 ^ ‘■’3 


W^-Wj WjW^ 


(25) 




If one choses the night time minimum density point in space 
defined by 


a 

V 

<5 

V 


Y + TT 
- 6 

s 


then the diurnal term B is equal to zero. Thus if three 
altitudes are chosen, h^, h^, and h^ , then the coordinates 

for the three altitudes where the diurnal term is zero are 
given hy 



= (R.-h^ 

) 

ans 

(Y + TT) 008 



= (H^.h^ 

) 

sin 

(Y + ’H') 000 



= (R.-h^ 

) 

sin 

(-6^) , 



( 26 ) 
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where R is the radius of the earth. Given the coordinates 
e 

for the .1 different altitudes, the sola- flux intensity, 
the geomagnetic index and the day of the year, one may deter- 
mine three densities from the Jacchia model. One then defines 
the function f^^ as 

f^ « (In i - 1, 2, 3 (27) 

where is the density determined from Jacchia evaluated 

with the ith position and xs found from evaluating equa- 

tion (23) with the ith position. If one then defines 

w^ - h^ for i = 1 , 2 , 3 (28) 

then coefficients of the night time vertical profile may be 
determined by equation (25) where 

a^ = for k = 1, 2, 3 . (29) 

The coefficients in the diurnal term may now be deter- 
mined by choosing the coordinates of the vehicle such that 
m 

^ is equal to 1 (daytime maximum density point). 

2 

This is the case when 



= (R.-^h^) 

C08 

Y (ioa 

6 

s 


^1 

= (R,-h^) 

ain 

Y aoa 

6 

s 

(30 ) 


= (R.-^h.) 

ain 

6 

s 




One then may define a new f^^ 

f^ = 1.1 (Pj/Pq) - (A^ + C^) for i = 1. 2, 3 (31) 
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where is the density determined from Jacchia evaluated 

with the ith position, and and are determined from 

evaluating equation (23) and equation (18) respectively with 
the ith position. If one then defines 

w^ - h^ for i - 1, 2, 3 (32) 


the diurnal coefficients may be determined from equation 
(25) where 


b,^ - Cj^ for k = 1, 2, 3 . (33) 

The proper altitudes to choose may be defined by the 
region in which one w'ishes to use the model. In the nume/ical 
experiments in Section 7, the following heights were chosen: 
hj = 150, h 2 ® 300, h^ = 450 km. If the model is to be used 
over a large range of heights then one may layer the nig^t 
time profile in two sections. 

By careful study of the Jacchia model a value of <{) » 37® 
(equation (12)) and m = 2.75 (equation (19)) has been adopted. 

(equation (22)) has the value of 1.224997 kg/m*. However 
all these values may be refined to obtain closer agreement to 
the Jacchia .nodel. 


7.0 Numerical Experiments 

* 

To compare the new density '.'odei to the Jacchia model , 
several numerical experiments have been carried out. Instead 
of comparing directly the computed densities, predicted sat- 
ellite positions are compared from a numerical orbit computa- 
tion program (Reference 4). The density model is used to com- 
pute the drag forces on the satellite. Several different or- 
bits were chosen for the comparison. Thus the position dif- 
ference is a good indication of the global difference between 
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the density models. The orbits chosen are the same as in 
Section 2.2.5, Reference 3. Also the ballistic number is an 
average value for the shuttle, BN ■ 100 lb/(ft)* , and the 
coefficients of drag, , is set to 2.2. Table I displays 

the chosen orbits and Table II gives the results. Also in- 
cluded are the results of the USSR and AMDB model. The 
hybrid model will henceforth be called the Babb-Mueller (B-M) 
model. Gus Babb, of FPD, originated the AMDB model and this 
author has Included the additional terms from the USSR model. 
The results in Table II are the differences in predicted posi- 
tion using the Jacchia model as compared to using the other 
three models or neglecting drag completely. 

The coefficients used in the USSR model correspond to 
an Fjq ^ = 75. This happens to be the actual F^^ ^ which 
occured in 1975. Thus one would expect that the USSR and 
Jacchia models to agree well in the year 1975. For all these 
cases in Table II the USSR model does show close agreement 
with Jacchia, but the now model (B-M) shows better agreement 
with Jacchia than any of the other models, including the USSR 
model. Case E has an epoch of 1977 which corrr>sponds to a 
Fjq ^ = 110. The USSR model still used the coefficients of 

Fj^q ^ '75. Thus one expects to see a disagreancuit to Jacchia. 
In Table II, one sees that this is the case. However, the B-M 
model still shows a close agreement. 


♦ 

In the Jacchia model, the daily changes in F._ have been 
neglected and also the hourly changes in the geomagnetic 

index (a ) are dropped. 

P 
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TABLE I: ORBITS USED IN PREDICTION EXPERIMENTS 



A 

B 

C 

D 

E 

perigee (km) 

220 

300 

300 

220 

Same as 

apogee (km) 

380 

600 

600 

380 

A except 

eccentricity 

.012 

.022 

.022 

.012 

epoch 

period (min) 

90.5 

93.6 

93.6 

90.5 

is 

argument of perigee 

0 

0 

180® 

0 

12:00 

ascending node 

0 

0 

0 

0 

January 1, 

inclination 

30® 

30® 

30® 

90® 

1977 

epoch 

12:00 

Janu 

ary 1 , 

1975 



TABLE II: POSITION DEPENDENCE ON DENSITY MODEL 


Time of 
integration 
( days ) 

No Drag 

Position difl 
USSR 

'erence (km) 
AMDB ] 

1 

1 B-M 

0.5 

1 

8.9 

0.2 

0.3 

.005 

1.0 

36.1 

0.8 

1.5 

.018 


(a) Orbit A 


TABLE II: CONTINUED 


Time of 
integration 
(days) 


Position difference (km) 


No Drag 

USSR 

AMDB 

B-M 

0.5 

.83 

.03 

.07 

.004 

1.0 

3.16 

.07 

.28 

.015 


(b) Orbit B 
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TABLE II: CONTINUED 


Time of 
integration 
( days ) 

Position difference (km 


No Drajt 

USSR 

AMDS 

B-M 

0.5 

.5 

< .001 

.01 

.004 

1.0 

1.8 

.02 

.04 

.018 


(c) Orbit C 


TABLE II: CONTINUED 


Time of 
integration 
(-days) 

Position difference (km) 

No Drag 

USSR 

AMDB 

B-M 

o 

8.6 

.1 

.3 

.008 

1.0 

34.4 

.5 

1.5 

.027 


(d) Orbit D 


TABLE II: CONTINUED 


Time of 
integration 

Position difference (km) 

No Drag 

USSR 

AMDB 

B-M 

0.5 

13.9 

3.4 

0.4 

.071 

1.0 

56.0 

19.2 

1.9 

. 275 


(e) Orbit E 







8.0 Conclusions and Outlook 


From the re.sults of the numerical experiments it is con- 
cluded that the B-M model gives better agreement to the 
Jacchia than either the USSR or AMDB model. This demon- 
strates the validity of the chosen approach. Additional 
studies may be conducted to determine the range of altitudes 
over which the model is valid and whether a layering of the 
atmosphere model is necessary. Also, additional studies need 
to be made to refine the angular parameter ( 4 >) which defines 
the position of the atmospheric bulge and the power exponent 
(m) in the diurnal term. 
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